Tunable gauge potential for neutral and spinless particles in driven lattices 
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We present a universal method to create a tunable, artificial vector gauge potential for neutral 
particles trapped in an optical lattice. The necessary Peierls phase of the hopping parameters 
between neighboring lattice sites is generated by applying a suitable periodic inertial force such that 
the method does not rely on any internal structure of the particles. We experimentally demonstrate 
the realization of such artificial potentials, which generate ground state superfluids at arbitrary 
non-zero quasi-momentum. We furthermore investigate possible implementations of this scheme to 
create tuneable magnetic fluxes, going towards model systems for strong- field physics. 



First introduced in electromagnetism, gauge fields play 
a central role in the description of interactions in physics, 
from particle physics to condensed matter. Currently 
there is a large interest to introduce gauge fields into 
model systems in order to study fundamental aspects of 
physics pp. Especially the emulation of synthetic elec- 
tric and magnetic fields for of ultracold atomic systems 
is crucial in order to extend their proven quantum sim- 
ulation abilities further, e.g. to quantum Hall physics or 
topological insulators. In this context, the analogy be- 
tween inertial and Lorentz forces triggered the simulation 
of homogeneous artificial magnetic fields using rapidly ro- 
tating trapped ultracold gases [2 J. Recently, several pro- 
posals ([3[ H] and references therein) and experimental 
realizations focused on the simulation of a gauge vector 
potential A either in a bulk system [SJ |6] or in optical 
lattices pj [8]. The realized schemes exploit the Berry 
phase which arises when the atomic ground state is split 
in several space-dependent sublevels, as in the presence 
of an electromagnetic field. Hence they rely on the cou- 
pling between internal and external degrees of freedom 
induced by laser fields. 

Here we demonstrate the generation of artificial gauge 
potentials for neutral atoms in an optical lattice without 
any requirements on the specific internal structure. As 
the realized scheme only relies on the trapability of the 
particle, it can be very widely applied to many atomic 
systems, in principle also to molecules and other com- 
plex particles. It is particularly interesting for fermionic 
systems, where in a many-body state governed by Pauli 
principle, the use of internal degrees of freedom often lead 
to conflicts with the creation of gauge potentials. As an 
important additional benefit, the internal degrees of free- 
dom of the particles can be addressed independently, e.g. 
by real magnetic fields or microwave excitations. 
In general, the presence of a gauge vector potential mod- 
ifies the kinetic part of the Hamiltonian describing the 
system. In a lattice, an artificial field can then be sim- 



ulated by engineering a complex tunneling parameter 
J = \J\ • e %e , where 6 is the Peierls phase. 
The central approach here is to control this phase via 
a suitable forcing of the lattice potential, acting at the 
single-particle level. We describe the general scheme for 
the creation of an artificial gauge potential and experi- 
mentally demonstrate the realization of a tunable vec- 
tor gauge potential in ID. We analyze the dynamical 
processes leading to the relaxation towards a superfluid 
ground-state at any desired, finite quasi-momentum. In 
a 2D lattice the described forcing can be used to create 
artificial magnetic fluxes smoothly tunable between — it 
and 7r. This opens the possibility to emulate strong- field 
physics with cold atom systems, in a regime currently not 
accessible in condensed matter because of experimental 
limitations on the magnetic field strength. 
For the general scheme, we consider a system of ultracold 
atoms in a deep optical lattice described by the Hubbard- 
type tight-binding Hamiltonian: 



H(t) 



(ij) 



"y]vi(t)hi + H 



(1) 



Here a\ , a i and hi denote the creation, annihilation and 
number operator for a particle (boson or fermion) of mass 
m localized at site i at r^ Jij > quantifies the tun- 
neling between neighboring sites, and ^on-site comprises 
time-independent on-site terms describing, e.g. inter- 
actions or a trapping potential. The system is driven 
by fast, off-resonant time-periodic potential modulations 
Vi(i) = Vi(t+T) of zero time average (^)t — ^ Jo dtvi(t) 
with period T. The case of a homogeneous inertial force 
F = —mx, which is created by shaking the lattice along 
the periodic orbit x(t) in space (see Fig. [TJa)), is de- 
scribed by Vi = —ri • F. Under the conditions de- 
scribed in Ref. [9 , one can show that the driven system 
is to good approximation described by the effective time- 




(b) 



F(t) 



(d) 



E(k) 




"l 2 



-►t 



-2n/d 



6/d 2n/d 



-■-►k 



FIG. 1. Periodic driving of the ID lattice, (a) Time-periodic 
movement of the lattice in real space, (b) Resulting periodic 
inertial force of zero mean. T1/T2 is the asymmetry parame- 
ter, (c) Realized complex tunneling elements on a ID lattice, 
with spacing d. (d) Effective single-particle dispersion re- 
lation. A complex phase induces a shift of E e ^{k) towards 
/c m in = 9/d (solid line). 



independent Hamiltonian: 
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The important terms here are the complex tunneling pa- 
rameters \J?f\e Wi i = Jiji&^-^/^T (see Fig. [ljc)), 

where Xiif) = ~ Jt Q dt' ' Vi{t f ) + (f*dt' Vi(t')) T - For sinu- 
soidal forcing, such a dynamic modification of tunneling 
is restricted to 6^- = or n [T0HT3] . It has recently 
been observed in several experiments [T4HT7] . We will 
now show that by suitable driving, the Peierls phases Oij 
can in fact be smoothly tuned to any value (8 € (0, 2ir]). 
This allows one to engineer vector potentials that can 
give rise to both artificial electric and magnetic fields. 
Similar ideas have recently been presented for a Bose- 
Fermi mixture in a triangular lattice [T8] . 

The Peierls phase 6ij can be varied smoothly, even 
though (xi)t = 0. Indeed (e l( ^ Xi_x ^)T can be complex 
provided the forcing breaks two symmetries that are also 
known to prevent ratchet-type rectification in classical 
[T9] and quantum [20] lattice systems, namely: (a) re- 
flection symmetry for a suitable time r, i.e. Vi(t — r) = 
Vi(—t — r) and (b) shift (anti) symmetry, i.e. Vi(t — T/2) = 
—Vi(t). This requires forcing with more than one fre- 
quency. The forcing function used in our experiment (see 
Fig. [lib)) induces a monotoneously increasing Peierls 
phase Oij = arg((e ima; -( r j- r ^ T ) with the amplitude of 
the force. 

The Peierls substitution directly links this phase to a 
gauge vector potential A via 0^ = J Vt dr • A/H (inte- 
grated along a straight path), which allows to create ar- 
tificial electric or magnetic fields. The modification of the 



Peierls phase in time results in an artificial electric force 
Fe = — A. It will be shown at the end of this manuscript 
that our scheme allows the emulation of a vector potential 
that gives rise to a finite artificial magnetic flux through 
the elementary plaquettes of 2D lattices. 
To demonstrate the power of this scheme, we have emu- 
lated a vector potential for a Rubidium Bose-Einstein 
condensate (BEC) in a one-dimensional lattice. The 
trapped particles are accelerated via frequency modula- 
tion of one of the lattice beams, a technique which is ex- 
perimentally straightforward. The resulting inertial force 
is comprised of a train of sinusoidal pulses separated by 
periods of rest of periodicity T=I\+T2= 1 ms as depicted 
in Fig. [lib). The zero mean value of this force prevents 
the transfer of a net-acceleration to the lattice. This leads 
to a renormalized tunneling matrix element of the form: 
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FIG. 2. Creation of complex tunneling matrix elements, 
(a) Absolute value of the tunneling parameter obtained from 
equation pf for our experimental parameters (Ti + T2=lms 
and T1/T2—2.I). (b) The measured Peierls phases in a ID 
driven optical lattice for different values of the forcing ampli- 
tude K are depicted as circles. The dashed red curve corre- 
sponds to the theoretically expected values (equation (pi)). 
(c) Quasi-momentum distribution of the BEC after 27 ms 
time-of- flight for different values of K. The Peierls phase 
as a function of K is deduced from the observed shifts of the 
interference patterns. 



where T\jT^ is the asymmetry parameter; K is the forc- 
ing amplitude (see Appendix) and Jo(K) the zero-order 
Bessel function. From equation ([3|, the induced Peierls 
phase in the effective tunneling element can be calcu- 
lated. For T2/T — >> 0, one recovers the harmonic driving 
introduced in [13] which only allows to control magni- 
tude and sign of J eff . In the opposite limit of small Ti/T, 
the tunneling amplitude | J eff | is only marginally affected 
whereas the phase depends linearly on the forcing am- 
plitude K. In order to avoid unwanted excitations of the 
system induced by strong forcing, we have chosen an in- 
termediate value T1/T2 = 2.1 for the experiment realized 
here. Note however that for such an asymmetry parame- 
ter the Peierls phase depends non-linearly on K (see Fig. 

lb)). 

The thereby introduced gauge potential A = hO/d mani- 
fests itself via a shift of the dispersion relation by 0/d as 
depicted in Fig. [TJd) . This follows directly from the effec- 
tive Hamiltonian pi), whose eigenstates are Bloch waves 
with the dispersion relation E(k) = — 2| J eff | cos(kd — 6). 
The gauge potential therefore allows for the genera- 
tion of superfluid ground-states (with group velocity 
v g = dE/dk = 0) at finite and tunable quasi- momentum 
k = A/h. As will be detailed later on, we experi- 
mentally observe the relaxation of the condensate quasi- 
momentum towards the minimum of the effective disper- 
sion relation. Therefore, the imprinted Peierls phase can 
be directly read out from the quasi- momentum distribu- 
tion revealed in time-of-flight after a sudden switch off of 
the lattice and the external potential. 
As a central result Fig. WCo) shows the experimental 
data together with the theoretical predictions from equa- 
tion pj. After increasing the forcing amplitude slowly 
(in up to 120 ms) to the desired value, the correspond- 
ing quasi-momentum distribution was recorded. From 
the obtained time-of-flight images, examples of which are 
shown in Fig. [2|c), we extract the Peierls phases 9 (see 
Appendix). We observe an excellent agreement between 
experiment and theory, thus proving the controlled gen- 
eration of an arbitrary vector gauge potential encoded 
into the Peierls phase G [0, 2tt]. In addition, the experi- 
mental images demonstrate the large degree of coherence 
maintained in the atomic sample throughout the shaking 
process. As an additional feature, Fig. [2|a) shows that 
the Peierls phase allows now to invert the sign of the ef- 
fective tunneling element without crossing |J eff |=0 via 
the rotation in the complex plane. 

In the following we will discuss the details of the relax- 
ation of the system towards non-zero quasi-momenta su- 
perfluid states, allowing for the described direct measure- 
ment of the Peierls phase. Note that for an homogeneous 
and non-interacting system, the initial Bloch wave at 
ki = remains an eigenstate of the effective Hamiltonian. 
Thus no transfer to states with k =^ is expected after the 
shaking is turned on. However, since we are working with 
interacting bosons and an external harmonic confinement 
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FIG. 3. Time-resolved measurements. (a),(b) A quench 
to a Peierls phase of -7r/4 within 1 ms induces an excitation 
of the dipole mode, leading to oscillations of the quasimo- 
mentum as depicted in (a) and observed in the time-resolved 
measurement of the quasi-momentum up to 0.5s.(b). The 
plain curve is a fit to a sinusoid from which one obtains 
/ ex p = 3.6 ±0.4 Hz, in perfect agreement with the expected 
dipole frequency of the dressed condensate. The final value of 
the Peierls phase is depicted as a dashed line, (c)-(f) Shift of 
the quasi-momentum towards a final Peierls phase of 37r/2 for 
a slow increase of the forcing amplitude over 160 ms (c),(d) or 
a sudden quench within 1ms (e),(f). The left graphs depict 
the center of mass dynamics in /c-space in both cases. On the 
right, series of time-of-flight images are reported as a function 
of time. The time-steps are 10 ms and 1ms, respectively. 



more effects come into play. 

When the gauge potential is ramped up from to Af, 
the condensate acquires a non-zero group velocity, re- 
flecting the presence of an artificial electric force Fe = - 
A. This velocity induces a displacement of the con- 
densate's center of mass position x c in the harmonic 
potential of frequency / (see Appendix). The result- 
ing restoring force induces oscillations both in position 
and momentum-space (see Fig. [3^a)). In Fig. |3|b), 
we report a time-resolved measurement of the conden- 
sate quasi-momentum after a quench to a final Peierls 
phase of -7r/4. The oscillations around the final quasi- 
momentum result from an excitation of the dipole mode: 
the measured frequency of 3.6 ± 0.4 Hz perfectly matches 
the expected dressed condensate frequency ^m/m* f for 
particles having an effective mass m* in the lattice of 
10 ±1 E rec depth with a tunneling amplitude of 0.3 Jbare 
(/theo=3.5 ±0.5 Hz). The coupling to non-zero quasi- 
momenta results thus from the underlying harmonic trap- 



ping potential. 

In addition, this center of mass dynamics is subjected 
to several damping mechanisms induced by the trap an- 
harmonicity or the lattice discreteness, which leads to a 
coupling to other collective modes and therefore to the re- 
laxation of the BEC towards the new equilibrium state. 
Therefore the duration of the ramp from to Af has 
to be compared with the timescale of those relaxation 
mechanisms. In Fig. [3] we compare time-resolved mea- 
surements of the quasi-momentum distribution for a slow 
ramp (Fig. [3Fd)) of A to a final Peierls phase #=37r/2, 
with a sudden quench (Fig. [3^f )) . As the gauge field is 
slowly increased, the BEC follows the shift of the disper- 
sion relation minimum, as depicted in Fig. [3Vc) . For the 
quench, on the contrary, for which the shift of the dis- 
persion relation occurs within one ms, the system can- 
not follow and thus relaxes into the nearest minimum 
of the effective band structure (see Fig[3^e)). For the 
chosen value this minimum lies on the left with respect 
to the original k=0 peaks and we thus find the BEC at 
k = -7r/2d. This demonstrates clearly that in presence of 
these relaxation mechanisms, the forcing does not induce 
a net particle current in the lattice unlike for ratchets 
but allows the engineering of ground-state superfluids at 
arbitrary non-zero quasi-momenta. 

Our scheme to generate complex tunneling offers fasci- 
nating possibilities to emulate gauge fields in higher di- 
mensions. In the following, few examples based on lattice 
shaking and modulated superlattices will be described. 
An artificial magnetic field is characterized by a finite 
magnetic flux through an elementary plaquette P of the 
lattice, given by the sum <£p = 6ij + Ojk + ••• + #/« 
mod 27r E (— 7r,7r] taken around P. Homogeneous forc- 
ing, implemented by lattice shaking, cannot give rise to a 
magnetic flux through plaquettes having pairwise paral- 
lel edges like square or hexagonal. It can, however, lead 
to a non-zero flux </>a through a triangular plaquette; 
<M = 0(\F o \d) - 20(\F \d/2) for F the amplitude of the 
force oriented parallel to one edge. As depicted in Fig. 
[4^a), continuously tunable staggered fluxes are realized 
(0v = — 4>a for the inverted plaquette). 

Our method enables one to continuously control the 
degree of frustration from zero (</>a = 0) to maximum 
(0a = tt) in lattice geometries like triangular or Kagome 
[2T| [22] . This control can be used to adiabatically pre- 
pare and study exotic quantum phases that are expected 
to appear in the strongly frustrated limit. Increasing 
the flux in a triangular lattice, frustration is reflected in 
the fact that the dispersion relation develops two min- 
ima that eventually, for 7r-flux, become degenerate (see 
Fig. [4^b)). In a triangular lattice with 7r-flux sponta- 
neous symmetry-breaking has recently been observed in 
an experimental simulation of classical magnetism [T7] . 
With the presented scheme, one is now able to study the 
influence of a small symmetry-breaking perturbation and 
to observe the many-body dynamics initiated by prepar- 




FIG. 4. Creation of magnetic fluxes in triangular lattices, 
(a) Staggered magnetic fluxes obtained in a 2D triangular lat- 
tice, (b) Dispersion relation for different values of the mag- 
netic flux |$| through a plaquette. 



ing a system in the upper minimum. 
Tunable magnetic fluxes through e.g. square or hexag- 
onal plaquettes can also be obtained via inhomogeneous 
forcing. A simple way of doing this, is to create an oscil- 
lating super-lattice as will be described elsewhere. 
Beyond the engineering of the hopping element, another 
essential feature of our scheme is the possibility to ad- 
dress the atomic internal degrees of freedom indepen- 
dently. Here a spin-dependent hexagonal lattice [23[ con- 
stitutes a promising system. In such a geometry the inter- 
play of next-nearest neighbor coupling, complex tunnel- 
ing matrix elements and spin are particularly interesting 
in connection with Haldane's model for instance [241125]. 
The very general method presented in this manuscript 
should allow for large sets of new gauge fields schemes 
for various particle classes in optical lattices. 
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LATTICE SHAKING: 

The lattice is accelerated along the z axis by modu- 
lating the frequency of one of the lattice beams with an 
acousto-optical modulator. We shake the lattice time- 
periodically with Cj/2tt = 1 kHz, as the modulation fre- 
quency has to be fast with rsespect to the atomic motion 
but small compared to band excitation energies. 
The inertial force acting on the atoms in the lattice frame 
consists of a train of sinusoidal pulses and is defined over 
a period as: 



F(t) 



Fq sin (uj\t) 




for < t < Ti = 
for Ti < t < T 



2tt 

2tt 



(4) 



The forcing parameter K is defined as: 
F d md 2 v 



K = 



fajj\ 



h 



(5) 



where v is the amplitude of the frequency modulation, 
ranging here from to 30 kHz. 

Many different analytical functions may be chosen for the 
time-dependent inertial force. 




APPENDIX 



ID OPTICAL LATTICE: 



The ID optical lattice is generated by two inde- 
pendent, phase stabilized laser beams at A = 830nm, 
intersecting at an angle of 117° in the xy plane, which 
results in a lattice spacing of <i = 486nm along z. The 
radial dimensions are only weakly confined by an optical 
dipole trap of trap frequency f z = 18 ± 2 Hz and f x = 45- 
55 Hz. We load a Bose Einstein condensate of Rubidium 
87 atoms into the lattice, which results in an occupation 
of around 80 pancake-like lattice sites. The lattice depth 
is 10±l£ rec , where E rec = 2ft 2 7r 2 /A 2 ra = 3.33 kHz. This 
results in a bare tunneling rate of roughly 0.007 E rec . 
The BEC is in the weakly interacting regime where 
U/J <<1, with on-site interaction U. 



FIG. 5. (a) Absorption images taken after 27ms time-of- 
flight over two periods of forcing. The center-of-mass oscil- 
lates under the effect of the forcing. The quasi-momentum, 
on the contrary, remains at any time described by the effec- 
tive Hamiltonian. (b) Plot of the fitted momentum of the 
background (resembling in good approximation the Wannier 
envelope), denoted kw, as a function of time. 



In order to induce complex tunneling parameters, the 
function should break the relevant symmetries (a) and 
(b) described in the text. For this, the frequency spec- 
trum of the force has to contain more than one compo- 
nent. 

A major advantage of the function denned in Q is the 
control of the sharpness of the Peierls phase dependency 
in K and the minimal amplitude of the effective hopping 
rate via the value of the asymmetry ratio T1/T2, where 
for the experiment T = T\ + T2 = 1 ms was chosen. 



TIME-OF-FLIGHT MEASUREMENTS: 

To reveal the quasi- momentum distribution of the BEC 
in the ID lattice, the lattice and the dipole trap poten- 
tial are suddenly turned off. We measure the absolute 
position of the interference pattern after 27 ms time-of- 
flight, with the quasi-momentum distribution obtained 
without forcing as a reference. When taking those pic- 
tures in the laboratory frame, one observes the quasi- 



momentum distribution of the state evolving accord- 
ing to the time-independent effective Hamiltonian, en- 
veloped by the momentum distribution of the lowest- 
band Wannier- function that oscillates like mx. Fig. [5] 
shows those oscillations over two periods of the forcing, 
corresponding to a Peierls phase of 0.64 ir. 
For the data presented in the text, the pictures were al- 
ways taken at times within a period where x — and 
averaged over 3 to 7 experimental runs. 



